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Abstract. We apply field theoretical renormalization group (RG) methods to describe the 
Tomonaga-Luttinger model as an important test ground to deal with spin-charge separation 
effects in higher spatial dimensions. We calculate the RG equations for the renormalized 
forward couplings g2R and g^R up to two-loop order and demonstrate that they don't flow 
in the vicinities of the Fermi points (FPs). We solve the Callan-Symanzik equation for Gr a \R 
in the vicinities of the FPs. We calculate the related spectral function and the momentum 
distribution function at p — kF.po — w. We compute the renormalized one-particle irreducible 
function T^ R (p,po = 0;A) and show it carries important spin-charge separation effects in 
agreement with well known results. Finally we discuss the implementation of the RG scheme 
taking into account the important simplifications produced by the Ward identities. 



1. Introduction 

The discovery of high-temperature superconductors and other related oxide based materials 
have attracted a lot of attention regarding the unique properties of low dimensional solids. The 
simplest and more notorious representation of a one-dimensional (Id) conductor is the the so- 
called Tomonaga-Luttinger (TL) modelpQ. That model essentially describes Id electrons in the 
vicinity of their right and left Fermi points interacting with each other by means of forward 
scattering processes. Due to the special equivalence relations between bosons and fermions 
and left and right charge conservation laws the TL model can be solved exactly by either 
bosonization[2] or, making full use of the Ward identities, by quantum field theory methodsjH]. 

The one-particle Green's function for the TL model can also be calculated exactly by those 
methods. However such a direct solution is provided only in coordinate space. To obtain the 
associated spectral function which determines the photoemission and the inverse photoemission 
spectra one needs to evaluate an intricate double Fourier transform and perform numerical 
approximations [I] . 

Bosonization methods are hard to implement in higher dimensions [Sj. In addition, those 
Ward identities simplifications are inherent to Id and to the exact number conservation of right 
(+) and left (-) particles. As a result to deal with strongly correlated fermions in two or more 
spatial dimensions one generally resorts to other schemes. One of those alternative approaches 
is the renormalization group (RG) method. RG methods were used with success to predict the 
various instabilities of the Id electron gas as a function of the existing coupling parameters [H]. 
The RG flows of the renormalized couplings are also in agreement with the bosonization solutions 
in several related problems. Despite that the RG method is not frequently used to perform a 



direct calculation of the single-particle Green's function and its associated spectral functional |Hj. 
Since in the TL model there are only low-energy bosonic spin and charge collective modes it is 
not clear that the right physics will emerge out of approximate fermionic RG schemes. Besides, 
it is also fair to say that RG methods have not yet been proved capable of dealing, in full 
force, with spin-charge separation effects. Within the g-ology framework the role of the forward 
coupling g4 has not been fully explored in the RG context so far. One should therefore try to 
establish how to deal with spin-charge effects in well known grounds, such as the TL model, 
before implementing RG methods in more difficult problems such as spin-charge separated states 
in higher spatial dimensions. 

In this work we apply the field theoretical RG method to describe the TL model in detail. 
This work is not intended to review exhaustively all the different RG applications for the TLM. 
We refer the reader to several authoritative review articles for that purpose 10 . Rather we 
concentrate here uniquely on the field theoretical RG description of the TLM. Using this scheme 
we derive the flow equations for the forward couplings g2R and g^R up to two- loop order. In 
calculating the respective renormalized one-particle irreducible functions and we show 
that non-parquet vertex contributions are the only source of logarithmic divergences in our 
perturbative calculations. Taking into account the self-energy corrections calculated earlier we 
demonstrate that, in two-loop order, the non-flow condition for g2R and g^R is assured by the 
exact cancelation produced by the contributions originated by the anomalous dimension and by 
the counterterms added to our renormalized TL Lagrangian. Based on our perturbative RG 
results we construct and solve the Callan-Symanzik equation for the renormalized one-particle 
Green's function Gr at the Fermi points. We show that it correctly possesses a branch cut 
structure due to its non trivial anomalous dimension. Using Gr we derive the spectral function 
and the momentum distribution function at the Fermi points. Our results are in agreement with 
earlier bosonization workj^j as well as more recent RG estimates [0]. Using a momentum RG 

(2) — l' 

scale parameter we use perturbation theory to derive T R = G R at the Fermi energy po = 
and in the vicinity of p = ±kp. One important feature emerges naturally from our result. There 
are now two emergent characteristic velocities which we relate immediately to the spin-charge 

(2) 

separation effects. Moreover, the general form of the resulting is in qualitative agreement 
with the Dzyaloshinskii and Larkin's Ward identity solution(DL)|S]. Our estimates are entirely 
perturbative since at that stage we only consider contributions up to two-loop order. To go 

beyond perturbation theory we relate the effective two-particle propagators and D ++ , 

introduced earlier by DL, to and and show how to implement the regularization, up 
to infinite order, of all the divergences which are originated in the perturbation series for these 
two one-particle irreducible functions. 

Recently new versions of the functional RG containing both fermionic and auxiliary bosonic 
fields have been developed with success [TH 115]. In one of those [Tl]. the Ward identities are also 
used to truncate the vertex functions hierarchical equations and as a result they are able to 
reproduce the correct TL spectral function at the Fermi points for the case in which g^ = 0. 

2. The TL Lagrangian Model 

Following the so-called g-ology convention the renormalized TL fermionic lagrangian is given by 
Ltl = I] (1 + AZ)iP R(a)(r (p,t)(id t - v F (\p\ - k F ))^ R{a)a {p,t) 

p,a=±,cr 

~ (92R + Ag2R)i> R ( +)ai (P3,t)i> R{ _ )rT2 (pi +P2 -P3,t)i> R{ - )a2 (p2,t)ijR (+)ai (pi,t) 
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where the momentum p is restricted to the interval —A < |p| — &p < A in the vicinities of the 
right and left Fermi points and AZ = Z — 1. In Figure ^ we display the direct couplings g2R 
and <74R. For practical purposes we can also take the limit A — > oo. We assume beforehand 
that neither the Fermi velocity nor the Fermi momentum are renormalized by interactions. The 
quasiparticle weight Z is nullified by interactions at the Fermi points since there are no stable 
quasiparticles in the TL regime but we assume tactically that we don't know this a priori. In 
general Z is the multiplicative factor which relate the bare and the renormalized fermion fields 
to each other. The counterterms associated with the contributions from Ag2R and Ag^R are 
included to regularize the logarithmic divergences which arise in perturbation theory. 

The corresponding right and left non-interacting Green's functions of this model are simply 
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Figure 2. The non-interacting Green's functions. 



Using the appropriate Feynman rules for Ltl we can now proceed with the perturbative 
calculations of the one-particle irreducible functions T^, and T^. In general, following the 
renormalization theory, we can relate those functions at the Fermi points to observable physical 
quantities. 

Due to the fact that only forward scattering processes are accounted for in Ltl the right (+) 
and left (— ) lines never mix in a single loop. As a result in this way no divergence arises from 
non-interacting particle-hole loop diagrams. Taking into account the spin summation factors we 
find instead 
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These n( )'s reduce simply to —i/itvp at the Fermi points if we take the limits qo 
q^O. 
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Figure 3. The non-interacting particle-hole bubbles. 



Another important feature due to the neglect of backscattering processes is that the 
symmetrized sums of all diagrams containing closed loops with more than two fermion lines 
vanish exactly 0. Thanks to this important simplification there are only divergent diagrams in 

the perturbative series expansions for both and in two-loop order. Those divergences, 
which are precisely the object of attention of renormalization theory, are cancelled out exactly by 
the counterterms added to the Lagrangian model. However at this order of perturbation theory 
we are also obliged to take into account self-energy corrections. Therefore before discussing 
the RG flow of the renormalized couplings we apply the RG method for the calculation of the 
self-energy S jR{+) . 



3. One-Particle Irreducible Function T 2 r up to Two-Loop Order 

Computing diagrammatically the self-energy X,r(+) (p,Po) for the TL model up to two-loop order, 
we find: 




Figure 4. The self-energy T. R ^ + y. 



Evaluating these diagrams we obtain 

(p,Po) = (P>Po) + (p,Po) + (p,p ) , (6) 

where 
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where Ap = p — k F , g iR = gm/v F and = 2v F X. 

Having established what £m_u is we move on to calculate the related one-particle Green's 

function Gy.ir = (t^J , ^ where 
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r H( + ) (P> Po) = PO - (P - fei?) - (p, Po) (10) 
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«(+) 



To determine AZ (&) we define such that, at the RG scale po = cj, and at the Fermi 
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point Ap = 0, rV + \ (p = kp^po = cj;cj) = a;. Using our two-loop results we get 

2 (»H-(fH") <"> 

Notice that the quasiparticle weight Z (u>) satisfies the RG equation 

din Z (uj ) , . 

u-^^=7, (12) 

where 7 is the anomalous dimension given by 

In view of this Z (u>) scales simply with the RG parameter as 

and since 7 > for non-zero forward coupling Z (uj) is nullified at the FP when ui — > in 
agreement with the fact that there are no well defined quasiparticles in the TL state. This 
result was obtained earlier in Refs. |1U1 lllj . 

4. Flow Equations for the Renormalized Couplings 

Following the renormalization theory we can relate the physical couplings g2R and g^R to the- 
particle irreducible functions and at the Fermi points. In this way for the external 
momenta pi = P3 = k F , P2 = ~k F and pio = P20 = P30 = <^/2 - Owe define 



r 2R (Pi = k F,P2 = ~kF]Pl +P2 -P3 = ~k F ,P3 = k F ;p w = p 20 = p 3 o = uj/2;u) = -ig 2 R (uj) 
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Figure 5. The channel. 



and similarly, for p l = p 3 = p 2 = k F and pi = P20 = P30 = v/2, 



(Pi = ^F,P2 = &f;Pi +P2 -P3 = ^F,P3 = &f;Pio = P20 = P30 = uj/2;uj) = -igm [u) 

(16) 

Using standard diagrammatic analysis as shown in the Figure we have at the vicinity of 
the Fermi point 
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To be consistent with the definition given above in ean !16l we then set 
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We display the diagrams for T^Ifp in Figure H3 Following the same strategy as before we 



indicate in Figure [7| how to proceed with the calculation of the counterterm — iAg2n( 
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Figure 6. The g 2 R channel. 
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Figure 7. The counterterm —iAg 2R (Lu) 



by 



From RG theory we know that the bare and renormalized couplings are related to each other 



g iB = Z 2 (to) (g iR + Ag iR ) 



(19) 



Since the bare couplings are independent of the RG scale lu they satisfy automatically the 
RG conditions 

dg2B dg 2 B n , orA 
lu— — = lu— — = (20) 
aw aio 

The flow equations for the corresponding renormalized couplings follow from these conditions. 
We find immediately that 
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Using our perturbative results for 7, Ag2 R and Ag^ R we then get that these flow equations 
are exactly nullified for the TL liquid: 
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A rigorous mathematical proof of the vanishing of the beta functions in all orders of 
perturbation theory in the TLM is given in Ref. jl.Sj . 



5. RG Calculation of the One-Particle Green's Function at the Fermi Points 

As we saw in the previous section the TL liquid is characterized by the non-flow of the 
renormalized forward couplings. This simplifies the RG approach to this problem. In this 
scheme we relate, quite generically, the renormalized Green's function G R to its bare analogue 
G B by 

G( a ) R (p,Po;92R;94R;uj) = Z~ l (Lu)G( a)B {p,po;g 2 B\g4B) (23) 

Since Gb is independent of the RG scale lu we derive the Callan-Symanzik equation (CZE)^B] 
differentiating G^b with respect to lu. Taking into account that the gi R s don't flow , in the 
vicinity of the Fermi points, for p = ±/ci?, the CZE reduces to 



(24) 



Since Gt a \ R has ordinary dimension of (energy) 1 , on dimensional grounds we have that 



d d 

^duj +P °dp^ ) G{a)R ^ = ±/cf 'P° ;w ) = - G (a)R(P = ±k F ,p ;Lu) 
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We can use this to rewrite the CZE in the form 



Po^~~ + 1 g (o)r(p = ±kF,Po;u) = (26) 

Considering that the anomalous dimension 7 is constant, this equation can be easily 
integrated to give 

Gu) R ijp = ±k F ,p ;u) = - ( — -) (27) 

which has indeed a branch cut structure as opposed to the simple pole nature, typical of the 
Fermi liquid one-particle Green's function. Consequently, for po < 0, at the Fermi points, 

G(a)R(P = ±k F ,p ;u;) = ~~ [rn J (cos tt7 + z sin ^7) (28) 
It follows from this that the spectral function A^ R (p = ±k F ,po;uj) reduces to 

A [a)R (p = ±k F , P0 ; u) = 6 (-po) ^ (29) 

This result is in agreement with earlier RG estimates of Kopietz and coworkers [S] . To calculate 
next the momentum distribution function n{p = ±k F ) at the FPs, it suffices to integrate ^4( a )_R 
over an energy interval of width 2oj around po = 0. We find 

, , sin7T7 

n(p = ±k F ) = -—<■ (30) 

ZTTJ 

which reduces to n (±fcp) = 1/2 in the limit 7 — > 0. 

Let us now consider the alternative limit po = and Ap = A ~ where the parameter A is 
now our RG momentum scale which will be used to take the physical system towards the FPs. 
If we repeat our RG analysis for making use of such RG scale we find at the FPs 

4wH = _»M|«M lD (|) pi) 

Since the forward couplings don't flow in the TL we can determine the corresponding Z (A) 
taking into account that we must have as before 

9™ n a dA 9iR 9iR H 9 2 2 R M /o 9 X 
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Consequently, 



and again we have that 
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as A — > for 7 > 0, which is consistent with the fact that there are no quasiparticles at the FPs. 
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Invoking once again the RG condition between the bare and renormalized 's we write 
r[+ } )fl (p, Po = 0;A) = Z (A) rg? )B (p, Po = 0) (35) 
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In the weak coupling limit in a regime consistent with our perturbative two-loop results Z (A) 
reduces to 



Z(A) = l + 7 h,(A) + .. 
Combining this with our earlier results produces 
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Consequently in the vicinities of the FPs 
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The presence of two different velocities violates the single-pole character of the one-particle 
Green's function and it is directly related to the spin-charge separation which takes place in the 
TLMfEj. This spin-charge separation is a feature which emerges naturally from our perturbative 
RG analysis. Notice that up to order O^r), m the weak coupling limit, vp ^/u c u s and r|^ R 
reduces to 
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(41) 



which is in qualitative agreement with well-known results [18] . However, our estimates for u c and 
u s are only valid up to two-loop order. For this reason they cannot be identical to the velocities 
obtained from the solutions of Schwinger-Dyson equations which include contributions up to 
infinite order 0. 



6. Ward Identities 

As demonstrated by Dzyaloshinskii and Larking the use of Ward identities greatly simplifies 
the analysis of the TL model. Making use of these identities and taking into account the 
simplifications produced by the cancelation of the symmetrized sum of all diagrams containing 
closed loops with more than two fermion lines we are naturally led, following DL, to introduce 
two interactions propagators D ++ and These two propagators can be used to describe the 

effective two-particle interactions in the TLM. D ++ and D_| are finite at the FPs and do not 

need to undergo any regularization procedure. In contrast with that, as we discussed in section 
4, the one-particle irreducible functions and produce non-parquet vertex contributions 
which need to be regularized order by order in perturbation series. Therefore to relate the DL 
interaction propagators back to the corresponding one-particle irreducible functions we must 
include appropriate vertex functions. Considering these feature we can readily extend our 
perturbative results writing the exact one-particle irreducible functions and in the 
form 



r 2R (Pl 3 P2,Pl + P2 -P3,P3) 



r+(P3)m>P3 -pi)iD + ^ {pi -p 3 ) 
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and 
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(43) 



These r^'s are represented diagrammatically in Figure |S1 Thanks to the already mentioned 

simplification inherent to the TLM, the propagators and D ++ are determined exactly by 

solving the Schwinger-Dyson equations indicated in Figures l$l and ITU1 Solving these SDE's for 

D ++ and , diagrammatically, we arrive immediately at the DL results for the charge and 

spin velocities, namely, u c = +54r/tt) 2 — ('Q2r/ 7T ) 2 an< ^ u s = V F- 

Since the interaction propagators are independent of renormalization parameters they are RG 
invariants. In contrast the vertex functions T's produce logarithmic singularities which need to 
be cancelled out by appropriate counterterms. As we showed before, in the RG framework, we 
solve this problem either by relating the bare and renormalized vertices to each other through 
relations of the type 



rf = z 2 r* = (i + Az 2 )r* 



(44) 



with a = +,— , or by constructing the necessary counterterm directly, order by order in 
perturbation theory. Following this latter route, up to two-loop order the renormalized vertex 
function is determined by the diagrams shown in Figure (jllj) . 

If we define T,? such that T? 1 = 1 it then follows that the vertex counterterm diagram 

which appears in Figure 1)12(1 is determined exactly by this condition. 
This gives 
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Figure 8. The exact one-particle irreducible functions and T^. 




Figure 9. The exact interaction propagators iD^ and iD 



Ar ?L = -(lr)M")= AZ (45) 

It follows from this that the vertex renormalization constant Z 2 must be identical to the 
quasiparticle weight: 

Z 2 = Z (46) 

Such an equivalence between the vertex renormalization constant Z2 and the quasiparticle 
amplitude Z was also pointed out in Ref. In fact, since the Ward identity(WI) is 

preserved by the renormalization process, this result emerges naturally from the WI itself. As a 
consequence of this simplification, at the FPs, the renormalized one-particle irreducible functions 

and reduce simply to the interaction propagators iD^ and iD ++ respectively, which 

by being RG invariants, produce vanishing beta functions to all orders in perturbation theory 
in agreement with the rigorous analysis of Ref. ^3] ■ 

7. Conclusion 

We implement the field theoretical RG method in the presence of spin-charge separation effects 
in one spatial dimension. We take the Tomonaga-Luttinger model as our reference test ground 
in dealing with future higher dimensional problems. As is well known trhe TLM was solved 
by bosonization and by quantum field theoretical methods which take into explicit account the 
simplifications produced by the Ward identities. Unfortunately those simplifications are inherent 
to one dimension and contrary to the RG scheme both methods are difficult to implement in 
more general situations. 

Using our RG method we calculate the self-energy up to two- loop order. From it we find a 
non-zero anomalous dimension and a quasiparticle weight which is nullified at the Fermi points. 
We compute the RG equations for the renormalized coupling constants g 2 R and g^R, up to 
two-loop order, and demonstrate that they don't flow in the vicinities of the FPs. We derive 
the Callan-Symanzik equation for the one-particle Green's function at p = kp,po = uj. The 




Figure 10. The propagators iD ++ and iD |_. 




Figure 11. The Renormalized vertex function r+. 



CSE is easily integrated producing an expected branch cut structure in Gr a \R. Using this we 
derive the spectral function and the momentum distribution function at p = &p. We calculate 

(2) 

the one-particle irreducible function (p,Po = 0; A) taking into explicit account spin-charge 

separation effects in the weak coupling regime. Our results are in qualitative agreement with the 
other approaches. Finally we discuss the inclusion of the Ward identities in the RG scheme. We 
show that their present forms are preserved upon the renormalization and thanks to them the 
one-particle irreducible functions in the vicinities of the FPs reduce to interaction propagators 
which can be solved exactly by appropriate Schwinger-Dyson equations. To conclude we add 
that some of the material discussed in this paper have appeared before in the literature in one 
form or another. We tried to relate, whenever possible, our results to some of those works. By 
presenting them here in a field theoretical self-contained form we hope to bring new insights 
concerning the RG applications with both spin-charge separation effects and Ward identities in 
more general problems. 




Figure 12. The vertex function counterterm 
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